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ABSTRACT 

A significant amount of energy induced by road unevenness and vehicle roll and pitch motions is 
usually dissipated by conventional shock-absorbers. In this paper, a novel active multiobjective 
control design methodology is proposed which explicitly includes, besides the usual control objec¬ 
tives on ride comfort, road handling and suspension stroke, the amount of energy to be harvested 
as an additional, though conflicting, control objective and allows the designer to directly trade-off 
among them depending on the application. An electromechanical regenerative suspension system 
is considered where the viscous damper is replaced by a linear electrical motor which is actively gov¬ 
erned. It is shown that the proposed control law is able to achieve remarkable improvements on the 
amount of the harvested energy with respect to passive or semi-active control strategies while main¬ 
taining the other objectives at acceptable levels. Simulative studies undertaken via CarSim are also 
reported that confirm the potentiality and flexibility of the proposed control design strategy. 
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1. Introduction 

In the last decade, the development and the analysis 
of novel control design methods for active or semi¬ 
active suspension systems for automobiles has received 
a great attention by academic and industrial researchers 
and many solutions based on different perspectives and 
approaches are now available (e.g. see Sabzehgar, Mara- 
vandi, & Moallem, 2014; Sammier, Sename, & Dugard, 
2000 and references therein). 

A typical architecture of an active suspension system 
consists of a spring in parallel with a viscous damper and a 
force actuator. A related technology is that of regenerative 
suspension systems (Montazeri & Soleymani, 2010), which 
aim at converting the kinetic and potential energy of the 
suspensions into electrical energy via an electromechani¬ 
cal generator instead of dissipating it. If such systems did 
equip a hybrid or electrical vehicle, the harvested energy 
could be used to recharge the vehicles battery and hence 
extend its mileage. In traditional vehicles, the harvested 
energy could be used to power the electrical subsystems 
and hence reducing the workload of the internal combus¬ 
tion engine. 

Most of the works on regenerative suspension sys¬ 
tems (see e.g. Montazeri & Soleymani, 2010; Pires, Smith, 
Houghton, & McMahon, 2013; Sabzehgar et al., 2014; 
Sagiv & Bobrovsky, 2011; Scruggs, 2010; Spelta et al., 
2010; Zuo, Scully, Shestani, & Zhou, 2010) have been 
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concerned with the technological and constructive details 
of the electrical motors and drivers and the control design 
aspects have found less consideration with rare excep¬ 
tions (see e.g. Scruggs, 2010). In fact, in the literature, the 
most used control strategy in such a kind of applications is 
the one which makes the electrical motor behave as a vis¬ 
cous passive damper. A more recent contribution in this 
respect can be found in Zhang, Cao, and Yu (2012). 

The aim of this paper is to propose a novel active mul¬ 
tiobjective Hoo control design methodology for regenera¬ 
tive vehicle suspension systems which explicitly includes 
the amount of harvestable energy among the usual 
control objectives on ride comfort, suspension stroke and 
road handling to trade-off. In the case of regenerative sus¬ 
pension systems, the presence of conflicting objectives, 
some to be minimised and some maximised, makes the 
design problem more difficult, giving rise to possible non- 
convex optimisation problems when approached directly 
within the Linear Matrix Inequalities (LMI) optimisation 
framework. 

This difficulty is addressed and solved in this paper by 
analysing and exploiting the relationship between elec¬ 
trical and mechanical power amounts in the regenera¬ 
tive system, flowing from the road to the electromechan¬ 
ical actuator. It will be shown that the employment of 
active control strategies in regenerative suspension sys¬ 
tems might improve considerably the harvested energy 
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Figure 1. Passive set-up (left)and regenerative set-up (right). 

w.r.t. passive or clipped strategies while satisfying, at the 
same time, the other control specifications. 

This paper extends preliminary results given in Di 
Iorio and Casavola (2012) and Casavola and Di Iorio 
(2012), where the MIPC control was introduced for the 
first time. Here, an exhaustive discussion and compar¬ 
isons with other more traditional control strategies are 
reported. 

7.7 The 2-DOF quarter-car model 

The 2-DOF quarter-car model depicted in Figure 1 is 
used for control design purposes. In the regenerative 
set-up, the passive viscous damper pertaining to the 
traditional suspension systems is replaced by an elec¬ 
tromechanical actuator that generates a force//) to suit¬ 
ably dampen the vertical motion. A classical quarter-car 
dynamical model for the passive set-up is given by 

m s z s (t) = -k s (z s (t ) - z u (t )) - c(z s (t) - z u (t )) ( 1 ) 

^u’^ui 0 — ks(z s (t) z u (t )) k u (z u (t) z r it)) 

+ c(z s (t) - z u (t )) ( 2 ) 

while the same equations can be rewritten as 


Note that, being the damper a dissipative device, it has 
been removed in the regenerative set-up considered here 
because the aim is that of maximising the harvestable 
energy. A force generator is used in its place to dampen 
the vertical motion and transform such energy in electri¬ 
cal power. 

7.2 Electromechanical actuator model 

In this paper, for simplicity, we consider a linear 
permanent-magnets electromechanical motor as force 
actuator (LPMA). The LPMA solution is preferable 
because the employment of magnets permits to avoid the 
excitation circuit that needs power for it to work. An alter¬ 
native solution may be that of using a rotary electrical 
motor with a mechanical transmission for motions con¬ 
version. However, this kind of solution has usually a lower 
efficiency because of the high frictions losses. The lin¬ 
ear dynamic equations of a generic LPMA actuator (Jang, 
Choi, Cho, & Lee, 2005) are given by 

di(t) 

v ai 0 = RaKt) + -F K v v(t ) (5) 

/ (0 — K v i(t) (6) 

where / is the damping force, i the armature current, v a 
the armature voltage and v the compression rate (verti¬ 
cal speed of the suspension stroke). The parameters R a 
and L a represent, respectively the armature resistance and 
inductance while K v is the force constant. It is very impor¬ 
tant to choose the parameters in such a way that the elec¬ 
trical power flowing through the machine is characterised 
by high voltages and low currents in order to have lit¬ 
tle Joule losses and, in turn, a good efficiency. A way 
to achieve this is to choose, compatibly to other physi¬ 
cal realisation constraints, actuators with high values for 
the force constant K v and low values for R a . Other kinds 
of losses are less important and are here neglected for 
simplicity. 


m s z s (t) = —k s (z s (t) - z u (t )) + f(t ) (3) 

Wl u Zuit) — k s (z s (t) z u (t )) 

- k u (z u (t) - z r (t)) - fit) (4) 

when the regenerative set-up is considered. In the above 
equations, m s and m u denote, respectively, the sprung and 
unsprung masses, z s and z u their positions and z r the road 
profile. In the passive set-up, the term c denotes the vis¬ 
cous coefficient of the damper and/(f) the force provided 
by the electromagnetic actuator while, in both set-ups, k s 
and k u denote respectively the suspension and the tyre 
stiffness. Finally, the quantity z s — z u represents the sus¬ 
pension stroke and — z r the tire deflection. 


7.3 Overall control architecture 

In Figure 2, the overall control architecture is presented 
which consists of two nested control loops ( Cascade Con¬ 
trol ), governed respectively by the Current Controller 
(CC) and the Regenerative Vibration Controller (RVC). 
CC is the internal controller of the four-quadrant DC- 
DC converter aimed at regulating the armature current 
loop so as to establish the correct bidirectional electri¬ 
cal power flows between the energy storage (battery) and 
the actuator. The circuital representation of the DC-DC 
converter is depicted in Figure 3. It is characterised by the 
following quantities: 
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Road profile 



Figure 2. Overall control architecture. 
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Figure 3. DC-DC converter circuit. 


SECOND QUADRANT 
power transfer from load 
to source 


FIRST QUADRANT 
power transfer from source 
to load 


The CC block consists of a PI control law which is 
used to maintain the armature current at suitable refer¬ 
ence values provided by the outer RVC controller. The 
related control signals conveyed to the DC-DC converter 
are modulated via a PWM. Its design however is not con¬ 
sidered here for brevity. Moreover, being its time con¬ 
stants expected to be remarkably smaller than that of 
the RVC loop, the inner current loop dynamics can be 
neglected during the RVC design. Then, hereafter we 
focus on the control design problem underlying the syn¬ 
thesis of the RVC controller which is expected to provide 
the usual control performance of standard active suspen¬ 
sion systems plus the ability to harvest as much energy as 
possible. 




THIRD QUADRANT 
power transfer from source 
to load 


FOURTH QUADRANT 
power transfer from load 
to source 


Figure 4. Four-quadrant characteristic for the DC-DC converter. 


• Vd : input dc voltage; 

• id : input dc current; 

• v 0 : output voltage whose amplitude and sign can be 
controlled; 

• i 0 : output current; 

• Ta +, Tb +, Ta Tb controlled switches; 

• Da +, Dg +, Da D# _: diods; 

• load : electromagnetic motor. 


2. Control design specifications and power flow 
analysis 

In classical active and semi-active control design prob¬ 
lems, three main design specifications are usually con¬ 
sidered: ride quality , road handling and suspension 
stroke. 

The ride quality is fundamental requirement for pas¬ 
sengers comfort and is mainly related to the sprung mass 
acceleration z s . A good suspension system has to main¬ 
tain this quantity within suitable limits at those frequen¬ 
cies where the human body is more sensitive to mechani¬ 
cal vibrations. For comfort evaluation the ISO-2631 (ISO, 
1997) standard maybe considered, where an index is pro¬ 
posed that gives a measure of passengers perception of 
vibrations: the ride index (RI) 


and behaves according to the four quadrants chart 

depicted in Figure 4. ride index := Ja 2 x + a 2 y + a 2 z , 


( 7 ) 
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with 


v^/ a2wi ^ dt i = x ^' z (8) 


a i(T) = J ~ 


P s (z r ,z u - z r ) 


Pm (f,Z s - Z U ) 


Pe(v,i) 



Figure 5. Power flows. 


where a w j is the acceleration along the ith axis weighed 
by the following filter: 

81.89s 3 + 796.6s 2 2 + 1937s + 0.14 

W (s) = ------ (9) 

s 4 + 80s 3 + 2264s 2 + 7172s + 21196 

The norm provides the above filter because it weights the 
frequencies whom the human body is particularly sen¬ 
sitive to. Note that lower values of the ride index imply 
better vibrations’ attenuation. In the present application, 
because only the vertical acceleration is of interest, the 
ride index essentially coincides with a z . Road handling is 
an another important objective and it is mainly related 
to vehicle and passengers safety. It becomes extremely 
important in sportive or race cars. It is usually represented 
by the tyre deflection — z r while the suspension stroke is 
usually represented by z s — z u . Obviously, suitable norms 
of the above control objectives are to be minimised for 
optimising the corresponding performance. 

A regenerative suspension system has of course to take 
into account all the above specifications for its design. 
However, an additional requirement on the amount of 
electrical harvested energy has to be considered too. To 
this end, relevant variables are the electrical power P e and 
the mechanical power P m of the electrical motor, defined 
as 


analysis of the regenerative suspension system is accom¬ 
plished in this subsection. From a systemic viewpoint, the 
latter can be considered as a series of two-port subsys¬ 
tems, as depicted in Figure 5 - with the dual variables on 
each terminal indicated between parenthesis - where P s , 
acting at the input terminal of the first subsystem, repre¬ 
sents the power exchanged with the road: 

P s (t ) := - k u z r (t)(z u (t ) - z r (t)) (13) 

In all terminals of the two-port system we adopt the con¬ 
vention that the power is entering into the system when 
positive. 

With reference to Figure 1, it is easy to understand that 
the road profile z r is the only exogenous signal in the sys¬ 
tem that excites the suspension and it is the only source of 
power to be harvested. In the forthcoming analysis, vehi¬ 
cles pitch and roll motions will not be considered for sim¬ 
plicity. Then, in the above simplified driving scenario no 
energy can be harvested on flat roads ( z r = 0). In order 
to exactly quantify the power flows we have to consider 
the total mechanical energy stored into the system, given 
by the sum of the kinetic energy of the masses and the 
potential energy of the springs: 


Pe(t) = rjP m (t ) (10) 

Pm(t) = /(f) (z s (t) - Zu(t)) (11) 

when the motor acts as a generator (that is P m (t ) < 0), 
P e (t ) is the useful power that can potentially be harvested. 
Note that P e (t) is a fraction of P m , with r] the correspond¬ 
ing conversion efficiency. On the contrary, when P m (t) > 
0, the electromagnetic actuator acts as a motor and the 
electrical power P e (t) is spent to produce the mechani¬ 
cal power P m (t). Finally, the power dissipated by the Joule 
effect can be expressed as 

P d (t) = RJ 2 (t ) (12) 

and motives the need of using high voltages and low cur¬ 
rents in this kind of applications. 

2.1 Power flows analysis 

In order to arrive to consistent guidelines on how to 
maximise the harvestable electrical energy, a power flows 


E(t) = 1 m s i 2 s (t ) + ^m u z 2 u (t) + 1 k s (z s (t ) - z u (t )) 2 
+ h«(z«(0 - z r (t)) 2 (14) 

Because we do not make use of dissipative devices like 
viscous shock-absorbers, the only way the power is 
exchanged and maybe dissipated is through the road and 
actuator terminals. The road irregularities can introduce 
energy or can absorb the potential energy stored in the 
springs when they discharge. Similarly, the actuator can 
absorb part of the energy stored in the system when it 
acts as a generator or can introduce energy into the sys¬ 
tem taking it from the battery when acts as a motor. 

Then, we can write the following instantaneous power 
balance: 

E{t) = P s (t) +P m (0 (15) 

where E(t) represents the instantaneous power stored 
into the system. By expanding expression (14) it is 
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found 

E(t) = m s z s (t)z s (t ) + k s (z s (t) - z u (t))(z s (t) - z u (t )) 

+ (O^u (0 

+ K(z u (t) - z r (t))(z u (t) ~ z r {t)) (16) 

and, by using the system dynamics (1) and (2), after direct 
manipulations one arrives to 

E(t) = f(t)(z s (t) - z u (t )) - k u z r (t)(z u (t) - z r {t)) 

(17) 

By direct examination, it follows that the first-term is 
related to the actuator whereas the second is related to the 
road. The above expressions related to P m and P s allow 
one to draw out several considerations: 

(1) An ideal control design would consist of having 
both quantities sign-defined for all time. In par¬ 
ticular, one would like 

P m (t ) << 0, P s {t ) >> 0 and P m (t) = -P s (t ), 

(18) 

conditions (18) would ensure that P s is large and 
always flowing from the road to the suspension 
system whereas P m is large as well and always 
flowing from the suspension to the battery. This 
means that all the power induced by road irregu¬ 
larities is harvested. In addition, condition P m (t) = 
—P s (t) would imply also that E(t ) = 0 which 
would ensure that the system has a perfect trans¬ 
fer of energy. However, in practical situations, only 
approximations of the above favourable control 
scenario are realistic. 

(2) Because f(t) is a directly manipulable variable, 
it is quite easy to make P m (t ) < 0 sign-defined. 
However, this choice is usually not convenient 
for harvesting a remarkable amount of energy 
because it does not take into account the role of 
P s into the play. In fact, in the favourable sit¬ 
uations whereby P m (t) < 0 and P s (t ) > 0, one 
arrives easily to the conclusion that the following 
condition 

[ \Pm(t)\dt < [ P s (t)dt (19) 
Jo Jo 

cannot be violated. In fact, Equation (19) results 
by taking the integral on any interval [0, T] of both 
sides of Equation (15) and considering that E(t) > 
0 for all t (£(0) = 0). Thus, it establishes a bound 
on the amount of electrical energy P e = r\P m that 
can be harvested and suggests that it might make 
nonsense to synthesise control strategies that try 


to maximise P m if P s is small and there is no way 
to increase its amount any further. 

(3) As far as P s is concerned, expression (13) rises 
several considerations. First of all, observe that 
P s is zero either when the road is completely 
flat ( z r = 0) or when the tyre deflection is zero, 
( z u (t ) — z r (t) = 0). While the first condition on the 
road flatness is quite intuitive the second is some¬ 
what unexpected and puts a trade-off between 
the harvestable energy and road handling perfor¬ 
mance. A further consideration is that P s depends 
on the regulated variable z u . Thus, if P s were made 
sign-defined and as large as possible, viz. P s > >0, 
the harvested energy could be increased via a suit¬ 
able control law regardless of the behaviour (the 
sign) of P m , that in this case need not be sign- 
defined (any active control law could be used). In 
fact, P s > 0 would provide a barrier for the internal 
potential energy of the suspension from flowing 
back towards the road. Then, even if some amount 
of energy would be provided by the actuator in 
some circumstances, this amount never would be 
lost because the system would result to be con¬ 
servative. These considerations are at the basis 
of the maximum induced power control (MIPC) 
design approach that will be presented in the 
next section and provides the main result of the 
paper. 

3. The maximum induced power control (MIPC) 
approach 

Based on the above considerations, it could be conve¬ 
nient to optimise P s by making it sign-defined and as large 
as desired. A rather simple heuristic will be presented 
hereafter which provides the basic idea of the approach. 
Consider expression (13). Then, P s could be made sign- 
defined ( P s > 0) if did exist a suitable control law able to 
enforce the following model reference condition 

z u (t) - z r (t) = -az r (t) (20) 

in closed-loop, where a is a scalar design parameter. If 
Equation (20) were true at each time instant the induced 
power would become 

P s (t) = k u az 2 r (t) (21) 

Therefore, the power coming from the road always would 
flow into the system and could be made as large as 
required by choosing suitable values for a. Obviously, 
it is not possible to make a arbitrarily large for having 
an arbitrarily large harvest energy. In fact, this would 
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increase the current i(t) and over a certain threshold the 
Joule effect would reduce the harvested energy instead 
of increasing it. In that case, P s would increase at the 
expenses of the energy taken from the battery The 
designer has to tune this parameter by increasing it until 
a convenient amount of harvested energy is achieved. 

Condition (20) provides a control design objective 
for energy harvesting maximisation. However, when the 
other, possibly conflicting, usual control objectives have 
to be considered the next problem is how to recast all 
of them in the same control design framework. To this 
end, because the T~L 00 control design approach has been 
largely considered in active and semi-active suspension 
control design problems we focus hereafter with this 
control strategy and it will be shown how the energy 
harvesting objective can be converted in the Tfoo-norm 
minimisation of a suitable causal transfer function. To 
this end, let us introduce the following performance 
output Zp 

z p (t ) = z u (t) - z r (t ) + az r (t) (22) 

and observe that z p (t) = 0 implies that the condition 
(20) holds true. Then, if a causal closed-loop transfer 
function between z r and z p did exist, say (with a little 
abuse of notation) z p = W Zp , Zr (s)z r , the energy/power of 
Zp could be minimised by minimising the Tfoo-norm of 
W Zp , Zr (s) provided that z r is in a specific class of signals 
with finite energy/power. Note also that, because both z r 
and z r appear in Equation (22), a causal transfer function 
W between z r and z p does not exist in general. On the 
contrary, it always exists when z r is used as the exoge¬ 
nous signal instead of z r , that is z p = W Zp ^ r (s)z r . Thus, 
in the next section a suitable state-space formulation will 
be described and the overall multiobjective T~L oo control 
design problem fully detailed. 

Finally, it is convenient to consider a as a stable filter 
rather than a scalar: this allows one to have system rep¬ 
resentations not in descriptor form, as they on the con¬ 
trary would result by the use of a pure constant Kf , and to 
have the energy harvesting prescription expressible as lin¬ 
ear combination of states, as required by the LMI control 
design method used here: 


constant r. In order the filter behaves practically as a con¬ 
stant, the time constant r should be chosen of course 
as small as possible without leading to numerical prob¬ 
lems and Kf as large as possible, to improve the energy 
harvesting performance. Higher order filters might help 
in achieving better performance but they complicate the 
control design. The analysis with higher order filters will 
be undertaken in future studies. 

3.1 Multiobjective / H OQ control design problem 

For causally realising z p in Equation (23) is necessary to 
work with a system description which takes into account 
as exogenous signal the derivative z r instead of z r . To this 
end, the state-space representation used in Chen and Guo 
(2005) is considered enlarged to include z r :=a(s)z r , 
with oi (s) the first-order lag (24). Then, the following 
state-space representation results 

x(t) = Ax(t ) + Bu{t) + Ez r (t ) (25) 

with state vector 

x [ z$ z u z s z u z r z u z r J (26) 

and system matrices 



0 1 

0 -1 

0 


k s 




-- 0 

0 0 

0 


m s 



A = 

0 0 

0 1 

0 


k s 

ku 



— 0 

-- 0 

0 


m u 

m u 





1 


0 0 

0 0 

— 




r _ 


0 

K v 


" 0 " 




0 


m s 



B = 

0 

, E = 

— 1 


K v 


0 


m u 


Kf 


0 


- r - 


Interestingly enough, the energy harvesting objective 
Zp(t) can now be rewritten as 


z p = z u — z r + a(s)z r 


(23) z p (t ) := z u (t) - z r (t ) + a(s)z r = x 3 (t) + x 5 (t) (28) 


A simple first-order filter 

a(s) = 


K 


f 


ST + 1 


(24) 


will be considered hereafter for simplicity with only two 
tuning parameters to be chosen: the gain Kf and the time 


The armature current u(t) = i(t) will be used as control 
command because we can neglect the DC-DC converter 
dynamics. 

According to the T-L 00 control design philosophy, let 
z(t) denote the vector of objectives to be minimised. 
As indicated, e.g. in Sammier et al. (2000), one typically 
includes in z(t) all specifications about ride comfort, road 








INTERNATIONAL JOURNAL OF CONTROL @ 747 


holding and stroke minimisation. Here, for simplicity of 
exposition, we consider only the ride comfort as stan¬ 
dard specification and add the maximisation of energy 
harvested and Joule s losses minimisation as new require¬ 
ments. To this end, the following objective vector can be 
defined 

Z(t) = [ z s (t ) i(t) z u (t) -z r (t) +a(s)z r (t)] T (29) 
that can be written as 

z(t ) — Cx(f) T Z)u(f) T Fz r (t ) (30) 

with 


k s 


~K V ~ 

-- 0 0 0 0 



m s 

, D = 

m s 

0 0 0 0 0 

1 

0 0 10 1 


0 


(31) 

As also common in the optimal FLoq control design 
approach, it is convenient to filter the objectives. To this 
end, by denoting W u (s) a low-pass filter for the current, 
we consider 


W z (s) = diag{W(s), W u (s), 1} (32) 

as the overall filter for the objective vector z, where W(s) 
has been already defined in Equation b(9). Also, an input 
filter W d (s) is used to shape the energy of the road profile 
z r . Then, by denoting with (A z , B z , C z , D z ) and (A d , B d , 
C d , D d ) the minimal state-space realisations of W z (s) and, 
respectively, W d (s), the overall augmented plant can be 
achieved: 

| X g (t) = A g x g (t) + B g2 u(t) + B gl z r (t) 

| z(t) = Cg\X g (t) + D g i 2 u(t) + D g nZ r (t) 


where 


A 0 EC d 


" ED d - 

B Z C A z B z FC d 

’ — 

B z FD d 

1 

o 

o 

a- 


B d 


B g2 — 


B 

B Z D 

0 


C g i = [ D Z C C z D z FC d ] 


(34) 

(35) 


D gU — D z FD d , D g u — D Z D (36) 


Then, by assuming for simplicity the state measurable 
(or estimable by a suitable observer) a state-feedback 
controller 


(37) 



Figure 6. Static state-feedback MIPC control scheme. 


can be used. Then, the following closed-loop system 
results 


[ x g (t) — C A g + B g 2 K)x g (t) + B g \Z r (t ) 
i ~ . (38) 

[ z(t) = (C g i + Dgi 2 K)x g (t) + D g iiZ r (t) 

and the control design consists in the minimisation of the 
FLoo -norm of the following transfer matrix: 


T (s) — (C g i + D g i 2 K ) (si — A g — B g2 K ) + D g n 

(39) 

This can be done, once suitable values for Kf and r have 
been selected, e.g. by solving the following LMI program¬ 
ming problem: 


min y 

X,Y,y 

~A cl +Al B g i ( C gl X + D n Y) T 

* -yi D* n 

* * —yl 

X = X T > 0 


< 0 


y > 0 


(40) 


where A c \ := A g X + B g2 Y for free matrices X = X T > 0 
and Y. If solvable (note that the feasibility of the above 
optimisation problem depends also on the choice of Kf 
and r in Equation (24)), the state-feedback controller is 
given by K= YX~ l and the overall control design scheme 
is depicted in Figure 6. 


3.2 Output feedback implementation 

In order to simplify the MIPC implementation, a state 
observer can be employed to implement the controller. 
The suspension displacement (z s — z u ) is sufficient to 
estimate the quarter-car models states via an observer. 
Unfortunately, because the derivative of the road profile 
is also to be used in the MIPC algorithm and one would 
like to avoid the direct measurement of it, the suspension 


u(t) = Kx g (t ) 
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Figure 7. Dynamic output-feedback MIPC control scheme. 


displacement measure is not enough informative for its 
estimation. On the contrary, the latter can be estimated by 
using additional signals, such as displacements and accel¬ 
erations. Note also that often such variables are already 
sensed in current vehicles, with sensors that are quite 
cheap and light. 

In the control literature many techniques have been 
proposed for real-time road profile estimation, most of 
them based on sensing the suspension stroke displace¬ 
ment, the sprung mass acceleration and the sprung mass 
displacement (see e.g. Doumiati, Victorino, Charara, & 
Lechner, 2011; Imine, Delanne, & M’Sirdi, 2006). These 
techniques are usually based on the use of a suitable state 
observers that exploit a dynamic model of the suspen¬ 
sion where the road profile and its derivative are included 
in the state vector. The method adopted here for assess¬ 
ment purposes is borrowed from Scherer, Gahinet, and 
Chilali (1997) and relies on the design of an output- 
feedback dynamic controller it is quite simple, being 
based on the basic quarter-car model, and achieves good 
performance. The required sensors are the suspension 
stroke displacement and the sprung mass acceleration: 
the sprung mass displacement is also a state compo¬ 
nent of the model but it is estimated through a dou¬ 
ble integration of the acceleration data. A trapezoidal 
method is used in Doumiati et al. (2011) for the inte¬ 
gration and gives reliable results. The control architec¬ 
ture is depicted in Figure 7 where the plant P has the 
same structure as in Equation (33) augmented with sen¬ 
sor measurements y(t) := [z s (t) — z u (t), z s (t)]. Note the 
advantage, remarkable for the implementation and con¬ 
trol design points of view, that neither z r (t) nor z r (t) 
compare directly in y(t): 


where 


Cg2 — \Cy, 02x5] 5 Cy — 


1 0 0 0 0 

0-^000 

m s 


Dg21 — [02xl] , D g 22 — [02xl] 


Please note that differently from the previous section, 
the measurement of the derivative of the road profile is 
no more needed. The controller will take the following 
form: 


x K (t) = A k x K (t) + B k y(t) , . 

1 u(t) = C K x K (t) + D k y(t) ^ } 

with matrices Ak, Bk, Ck, Dk determined by solving the 
following LMI-based problem: 


min y 

X,Y,Ak,B k ,Ck,Dk,y 

R\ + R\ R 2 R 3 

★ —yl (D 2 i +D 2 2D^Dn) T 

★ ★ —yl 

X = X T > 0 
Y = Y t 


(43) 


where 


AY + B g 2 Ck A + D^Cy 
Ak XA + 


Ri 


B g i + Bg 2 DKDgii 

XB g2 + B K D g2 2 


R 3 — ^ C Z Y + Dg22^K C Z Y + D g 22D K C y 


If the previous optimisation program admits a solution, 
the matrices of the PL 00 controller (42) can be obtained 
as 


D k = D k 

C k = (C K ~ D K C y Y) (M t ) _1 
B k = N~ 1 ( 'b k -xb u d k ) 

A k = N~ l (A K - NB K C y Y - XB U C K M T 
-X(A + B u D K C y )Y)(M T y 1 

where matrices M and N are computed by exploiting the 
factorisation NM T = I — XY. 


x g (t ) = A g Xg(t) + B g2 u(t ) + B g iZ r (t) 
z(t) = CgiXg(t ) + D g i 2 u(t) + DgnZ r (t) (41) 

y(0 — Cg2 Xg (f) Dg22 u(t') Dg2\ Z r (t ) 


4. Simulation results 

Several simulations have been performed for testing the 
effectiveness of the proposed approach. Power spectral 
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Table 1. (a) Vehicle data and (b) actuator data, from 
‘Bautzs, Danaher Motion'. Model: MA508G*. 


Parameter 

Value 

(a) 

Sprung mass 

1274 Kg 

Front wheels'mass 

35.5 Kg 

Rear wheels'mass 

35.5 Kg 

Roll inertia 

606.1 Kg • m 2 

Pitch inertia 

1523 Kg • m 2 

Front suspension stiffness 

27 kN/m 

Rear suspension stiffness 

30 kN/m 

Front tyre stiffness 

228 kN/m 

Rear tyre stiffness 

228 kN/m 

(b) 

Weight 

16 Kg 

Length 

36.5(+20) cm 

Force constant 

938 N/A 

Inductance 

5mH 

Resistance 

0.05 Q 


*This is an electrical servo actuator model. Retrieved from: 
http://www.wiki-kollmorgen.eu/wiki/DanMoBilder/file/ 
pdf_archiv/bautz/linearactuators/ma_e.pdf 


density of the roads profiles complying with the ISO- 
8608 standard has been generated and used so as to make 
the analysis repeatable. The vehicle chosen for the exper¬ 
iments is a standard class-C car, whose parameters are 
listed in Table 1 together with the actuator s data. 

4 .1 MIPC control simulations 

A output-feedback dynamical MIPC control law was syn¬ 
thesised according to Section 3.2, in three different ver¬ 
sions, using three different tunings of the weighting fil¬ 
ters and parameters (Kf, r). In this respect, the YALMIP 
(Lofberg, 2004) software tool has been used for the com¬ 
plete design of the controllers. The first version (MIPC 
comfort) has been designed for maintaining the perfor¬ 
mance very close to those pertaining to a classical passive 
suspension system with minimal requirements in terms 
of harvested energy. This first controller corresponds to 
the choice (Kf = 13, r = 0.0005) where the value of r 
has been selected by means of an exhaustive trial-and- 
error tuning procedure. A second version (MIPC trade¬ 
off) has been synthesised with increased requirement in 
terms of harvested power but with acceptable degrada¬ 
tion on the other objectives. This second controller cor¬ 
responds to the choice (Kf = 15, r = 0.0005). Finally, a 
third version (MIPC power) has been designed by max¬ 
imising the harvested power with (Kf = 18, r = 0.0005). 
This last choice is finalised to drain into the system large 
amounts of power, more than those that the system can 
handle without severe performance degradation on the 
other control specifications. It has to be considered as a 
limit solution for this system in looking for active con¬ 
trol laws that maximise the amount of harvestable power. 
Note however that higher values for Kf could have been 


chosen in this example, allowing a further maximisation 
of the harvestable power. 

For comparison purposes, we have also implemented 
a regenerative damper control law (44), quite simple and 
common in the literature (Goldner, Zerigian, & Hull, 
2001), where the electromagnetic unit is governed in such 
a way to mimic the behaviour of a passive viscous damper, 
with the damping parameter c optimised for maximising 
the harvested energy, c = 2051.4 ^ in this example: 

i(t ) = — (z s (t) - z u (t)) (44) 

4.7.7 Frequency behaviour of MIPC 

It is very interesting to look at the frequency behaviour of 
the MIPC controller. Harvestable power is defined by two 
main variables: — z r and z r . When the phase between 
these variables is 180°, P s is positive. For these reasons, 
the transfer function to be analysed is the one between 
the derivative of the road profile and the tyre deflection: 

Z u (s) — Z r (s) 

-=- (45) 

Z r (s) 

which gives an indication of the harvested power con¬ 
tribution at the different frequencies. The Bode plots of 
the above transfer function for different values of the 
weighting terms have been depicted in Figure 8 for the 
RD and MIPC control strategies. In particular, it is possi¬ 
ble to observe how the MIPC tuned for trading-off tends 
to amplify z u — z r especially at that frequency where the 
phase is exactly 180°. Observe also how the response for 
the comfort-tuned MIPC is quite similar to that of the RD 
control strategy. 

An explanation is that MIPC tends to amplify the mag¬ 
nitude at all frequencies if the maximisation of the har¬ 
vestable power is the only considered control objective. 
On the contrary, when also the ride quality is taken into 
account in the control design, the power is mainly har¬ 
vested around frequencies where the phases is 180°. In 
fact, in such a cases only the resonance is amplified. 

4.2 Simulations on the quarter-car model 
4.2.1 Bump of road 

In the first simulation, the quarter-car model has been 
excited by a 15 x 150 mm bump at a speed of 50 ^ 
depicted in Figure 9. The system responses related to 
the action of the different controllers are depicted in 
Figures 10-12. In particular, in Figure 10 the relative 
suspension deflection is displayed in order to verify that 
it remains within reasonable suspension limits. The Root 
Mean Square (RMS) value of the vertical acceleration 
of the chassis a z (t) has been reported in order to asses 
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Bode Diagram 



Figure 8. Frequency response of Z “ ( ^ ^ {s) with MIPC. 
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Figure 9. Bump of road of 15 x150 mm at a speed of 50 

the comfort experience during the simulations. Details 
about the CC in terms of current reference i(t) sent to the 
current controller CC, control signal and current track¬ 
ing error, respectively, have been depicted in Figure 11. 
Finally, in Figure 12 the harvested power for each 



- Regenerative Damper 



i 


IV/IIDr' nnmfru't 

MIPC trade-off 


1 

K 







1 

r 



0 0.5 1 1.5 2 2.5 3 



Figure 10. Bump of road test: chassis displacement and vertical 
acceleration. 


controller has been reported. It is clear that the MIPC- 
based scheme is able to gain much more power with 
respect to the RD approach. 

4.2.2 Different road classes 

Figure 13 has been obtained by Simulink simulations on 
the quarter-car model (for a single suspension). It depicts 
in a log scale how the harvested average power changes 
with the changing of vehicle s longitudinal speed and road 
roughness. In the picture, three groups of four plots are 
reported. The three groups correspond to three different 
road classes (from the below: type A, type C and type 
D) according to the ISO-8608 standard. Each group of 
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Figure 11. Bump of road test: signals related to the CC. 
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Regenerative Damper 
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Figure 12. Bump of road test: electrical power under different controllers. 
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Figure 13. (Colour online) Average power for different controllers 
and road types versus longitudinal speed. 

four plots (RD - blue, MIPC comfort - red, MIPC trade¬ 
off - green, MIPC power - black for online version of 
the paper) reports the average harvested power resulting 
from driving at constant speed on a straight road under 
the four controllers. 

It is easy to see that there is a remarkable gain in the 
harvested power by using the proposed MIPC approach 
with respect to the use of the RD control, the amount 
of which depending on the tuning of the weighting fil¬ 
ters. For example, consider the upper-limit of the C class 
road for a speed of 50 km/h (typical urban situation): 


under RD, the average power is about 40 W, while for the 
three versions of MIPC are 50,80 and 160 W, respectively. 
Noticeable enough, this is the harvested power from a sin¬ 
gle suspension (single wheel). 

Finally, plots of the electrical powers gener¬ 
ated/consumed by the actuator are reported in Figure 14. 
Under the passive RD control, the electrical power is 
always negative as expected: this means that the actuator 
acts always as a generator and the power flows always 
from the road to the battery. On the contrary, the active 
behaviour of the MIPC control law is evident: the power 
is sometime absorbed and sometime used to correct 
dynamic behaviour of the suspension system. Further¬ 
more, both the instantaneous peak values (which are 
not always negative) and the average amount (which is 
always negative) of the electrical power are higher than 
in the RD case. 


4.3 Simulations on CarSim 8 environment 

More detailed simulations were undertaken in the Car¬ 
Sim 8 environment. The controllers were tested on a 
full (four wheels) vehicles model, with detailed nonlin¬ 
ear dynamics. Four regenerative suspension systems were 
used for the four wheels, each one controlled separately 
(in a decentralised way) by the controller under test. 
Investigations on the use of a multivariable centralised 
controller will be undertaken in the future studies. 
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Electrical Power with Regenerative Damper 
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Figure 14. Electrical power under different controllers. 


Table 2. Numerical results for tests I and II (four 
suspensions). 


Controller 

Rl 

Max tyre defl. 

Av. Elect. Power 

[Test 1] 

RD 

0.55 

2.09 cm 

190.7 W 

MIPC (comfort) 

0.65 

3.3 cm 

265.5 W 

MIPC (trade-off) 

0.89 

4.4 cm 

286.5 W 

MIPC (power) 

2.89 

8.5 cm 

391.6 W 

[Test II] 

RD 

0.9 

2.63 cm 

495.35 W 

MIPC (comfort) 

1.38 

3.5 cm 

550.6 W 

MIPC (trade-off) 

1.53 

6 cm 

740.8 W 

MIPC (power) 

2.89 

13.8 cm 

1233.5 W 


In the first test, the car drives on a straight road with 
average road roughness (ISO road class C) at 70 km/h 
while in the second it drives at 18 km/h on a poor road 
(ISO road class D). Table 2 reports the corresponding val¬ 
ues for comparisons and performance evaluation. 

In the third and last test, the car drives at 70 km/h on 
a twisty road with low roughness (ISO road class A). The 
aim of this simulation is to show how the roll and pitch 
motions contribute to power generation. This could be of 
interest even if their motions have not been optimised by 


Table 3. Numerical results for test III (four suspensions). 



RD 

MIPC 

(comfort) 

MIPC 

(trade-off) 

MIPC 

(power) 

Ride index 

0.125 

0.23 

0.25 

2.4 

Electr. Power [W] 

47.1 

64.3 

65.9 

72.6 

Max roll [deg] 

[-4,33] 

[-1.9,2.7] 

[-1.9,2.7] 

[-17,2.5] 

Max pitch [deg] 

[-07,2.2] 

[-0.02,1.6] 

[-0.02,1.7] 

[-0.02,1.4] 

Max tyre defl. [cm] 

1.2 

1.97 

2.1 

2.8 


the four decentralised controllers. In this respect, Table 3 
shows the achieved numerical results. 

For the sake of clarity, the videos of these simula¬ 
tions have been published on the web page £ http://t-con- 
casavola_diiorio_tedesco.kissr.com/ > . 

In order to analyse the chassis’ motions effect on the 
amount of harvestable energy, it is very useful to look at 
the plot of the energy versus time corresponding to the 
use of MIPC with comfort tuning in Figure 15. The con¬ 
stant road roughness’s power contribution is depicted by 
the green line having a constant slope (the slope repre¬ 
sents the instantaneous power). During the simulations, 
some interesting events occur: first, after about two sec¬ 
onds, the car starts cornering a high rolling bend. During 
this phase, the slope of the blue energy graph increases 
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Figure 15. Energy plots on test III (the instantaneous power is rep¬ 
resented by the slope of the plots). 


Table 4. Power contributions in test III (four 
suspensions). 


Source 

Power contribution [W] 

Road roughness 

14.96 W 

Roll (during the first turning) 

290.64 W 

Pitch (during the braking phase) 

76.33 W 


remarkably. Between time instants 7 s and 10 s a straight 
stretch occurs, where only the road roughness contributes 
to the energy recovery. Finally, around time instant 11s, 
the driver claps sharply on the brake and the pitch motion 
is excited. In this brief time interval, the harvested power 
increases too. 

The average harvested power during the entire simula¬ 
tion is 64.33 W. By evaluating the parts of it at those time 
intervals when a specific event occurs and by subtracting 
the constant road roughness contribution, it is possible to 
make an estimate of how the chassis' motions contribute 
to power generation (see Table 4). 

Remark 4.1: From the above analysis, a consistent 
amount of energy seems to come from roll and pitching 
motions. Note that the proposed design framework has 
not been optimised for maximising the energy harvested 
during such manoeuvres where the harvesting would be 
more profitable. For space limitation such a further analy¬ 
sis will be presented elsewhere. It is expected that it can be 
achieved by redesigning the control strategy on a full-car 
vehicle model. 

Remark 4.2: For comparisons about the ride comfort 
performance pertaining to the various control laws, 
please observe that the above tables show the ride index 
of the RD control strategy. This measure gives a good way 
of comparisons, because RD behaves like a typical passive 
system. Note also that the fact that the proposed approach 
generates less performance on ride comfort compared to 
passive suspension, even for the comfort case, is a design 


choice. However, the MIPC method is flexible enough to 
achieve better ride comfort than RD when required, e.g. 
in applications where the ride comfort is considered as 
the first priority and the energy harvesting is just a by¬ 
product to achieve if possible. See Casavola, Di Iorio, and 
Tedesco (2014) for some preliminary results on this class 
of strategies. 

5. Conclusions 

The usage of efficient and light regenerative suspension 
systems can contribute considerably in modern electri¬ 
cal/hybrid cars to their autonomy . The integration of 
such systems with other existing devices, as regenerative 
brakes, may help the diffusion of this kind of cars, provid¬ 
ing together a total amount of many tens/hundreds watts. 
This paper has presented an TLoo -based output feedback 
scheme for the control of a suspension system in the 
presence of an electromechanical actuator able to harvest 
mechanical energy. 

An important contribution of this work is the demon¬ 
stration of the existence of a trade-off between road 
handling and energy harvesting specifications that lim¬ 
its the design of every kind of regenerative suspension 
system. It has been also demonstrated that the perfor¬ 
mance/harvesting trade-off can be dealt with effectively 
by active control laws as the MIPC here proposed. By 
classical control techniques, it is possible to design a reg¬ 
ulator that achieves classical performance specifications, 
but it is difficult to act directly on the energy harvesting 
requirement. The proposed algorithm (MIPC) gives to 
the designer a way to handle the trade-off in an appro¬ 
priate manner, depending on the application. 

Future works will concern the use of multivariable 
MIPC control laws designed on the basis of a full four- 
track model in order to further optimise the performance. 
Other studies related to the enhancement of the flexi¬ 
bility of the scheme will concern the use form of gain¬ 
scheduling MIPC control laws, mostly achieved via fuzzy 
or LPV control design methods. 
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